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O ■ ABSTRACT 

o; 

We analyze the linear stability of a dilute, hot plasma, taking into account 
the effects of stratification and anisotropic thermal conduction. The work is mo- 
tivated by attempts to understand the dynamics of the intracluster medium in 
C*") \ galaxy clusters. We show that magnetic field configurations that nominally sta- 

f— ,| bilize either the heat-flux driven buoyancy instability (associated with a positive 

thermal gradient) or the magnetothermal instability (negative thermal gradient) 
can lead to previously unrecognized g-mode overstabilities. The driving source of 
the overstability is either radiative cooling (positive temperature gradient) or the 
JL \ heat flux itself (negative temperature gradient). While the implications of these 

overstabilities have yet to be explored, we speculate that the cold fronts observed 
in many relaxed galaxy clusters may be related to their non-linear evolution. 

Subject headings: instabilities — plasmas — magnetohydrodynamics (MHD) - 
00 ' convection — galaxies: clusters: intracluster medium 

in ■ 

^r ; 

1. Introduction 

o: 

The dynamical stability properties of a dilute plasma change dramatically when even a 
weak magnetic field with no "J x B" force is present. This is a consequence of anisotropic 
^ transport: if the mean free path exceeds the relevant gyroradius, heat and (angular) momen- 

tum are transported only along magnetic lines of force (Balbus 2000, 2004). In particular, 
in a stratified atmosphere of dilute plasma, the Schwarzschild stability criterion (i.e. specific 
entropy must increase upwards in a gravitational potential) no longer governs convective 
stability (Balbus 2001). Instead, the presence of any temperature gradient is potentially 
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convectively unstable, depending on the relative alignment of the magnetic field and ther- 
mal gradient (Quataert 2008, hereafter Q08). The only configurations that are nominally 
stable are i) an upwardly decreasing temperature profile whose gradient is precisely aligned 
with the magnetic field, and ii) a upwardly increasing temperature profile whose gradient 
is precisely orthogonal to the magnetic field. Any other configuration is unstable. In this 
Letter, we show that even these "stable" configurations are in fact subject to overstability. 
Interestingly, the non-linear evolution of the thermal gradient instabilities tends to drive the 
system towards these overstable states. In essence, a dilute stratified plasma that is hot at 
its core wants to let the heat escape; a cold-core system wants to self-insulate itself. The hot 
core instability is known as the magnetothermal instability (MTI). The cold core instability 
is known as the heat flux buoyancy instability (HBI), because the free energy source is the 
thermal flux itself. 

Both the MTI and HBI may be of considerable importance for understanding the be- 
havior of the intracluster medium (ICM) in clusters of galaxies. The outer regions of ICM 
atmospheres have a negative temperature gradient and MTI-driven turbulence may be im- 
portant for energy and metal transport (Parrish, Stone & Lemaster 2008). Furthermore, in 
clusters with cooling cores, radiative losses lead to the formation of a positive temperature 
gradient in the central regions of the ICM: a cool ICM core is embedded in a hotter ICM 
atmosphere (e.g., review by Peterson & Fabian 2006). A longstanding problem has been to 
understand the thermal stability of the cool inner regions of ICM atmospheres. In particular, 
why do radiative losses not cause run away cooling in these regions? The standard paradigm 
is that regulated energy injection into the ICM by a central active galactic nucleus (AGN) is 
a stabilizing influence. However, the finding of Q08 the positive temperature gradient within 
the ICM of cooling core clusters could be convectively unstable suggests that the destabilzing 
role played by thermal conduction merits closer scrutiny in this context. 

A major difficulty with invoking thermal conduction to stave off radiative core collapse 
of the ICM (e.g., Bertschinger & Meiksen 1986, Zakamska & Narayan 2003) has been the fine 
tuning required (Conroy & Ostriker 2008). A little too much thermal conduction evaporates 
the cool core; too little allows the collapse to proceed essentially unimpeded. Balbus & 
Reynolds (2008) hypothesized that the combination of thermal and radiative losses may be 
self-regulating. In this view, suppression of the thermal conduction is not complete. Instead, 
through field line tangling and a possible reverse HBI convective thermal flux, the radiative 
core is able to draw in the requiste heat flux for marginal stability. Global MHD models 
of cluster cores have recently tested these ideas (Bogdanovic et al. 2009; Parrish, Quataert 
& Sharma, 2009, 2010; Ruszkowski & Oh 2010). These simulations show that in idealized 
static atmospheres, while there is significant heat transport into the cooling core, after several 
cooling times the non-linear development of the HBI wraps the magnetic field onto shells, 
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insulating the core from further conduction and fostering thermal collapse. Models with 
sufficiently strong forced turbulence (perhaps the effects of an AGN, or stirring by galaxy 
motions and/or sub-cluster mergers) can prevent the field-line wrapping and stabilize the 
cluster. However, such models predict temperature decrements in the core of at most a factor 
of two; real cooling core clusters often have much steeper temperature gradients. Moreover, 
the forced turbulence is assumed to be volume-filling, the justification for which is unclear. 

In this Letter, we return to the fundamental question of the linear stability of an at- 
mosphere of dilute weakly magnetized plasma in the presence of a background temperature 
gradient. We generalize the treatment of Q08 to include optically-thin radiative losses. Our 
principal result is the finding of new overstabilities of dynamical waves. More precisely, 
we find that nominally stable configurations resulting from the non-linear evolution of the 
HBI (i.e., temperature increasing upwards and magnetic field essentially horizontal) generate 
overstable g-modes via radiative losses. Nominally stable configurations resulting from the 
non-linear evolution of the MTI (i.e., temperature decreasing upwards and magnetic field 
essentially vertical) generate overstable g-modes via anisotropic thermal conduction. In addi- 
tion to furnishing a more complete formal picture of the stability properties of dilute plasma 
atmospheres, these findings may have significant implications for the physical behavior of 
the ICM, and should guide future simulations. 

In the next section, we present the calculation in detail, and in the final section of this 
Letter, we conclude with a brief discussion of the implications of our findings. 



2. Analysis 

We use the standard equations of MHD, with the entropy equation agumented with 
anisotropic thermal conduction along magnetic field lines (Braginskii 1965) and radiative 
losses (e.g. Field 1965). The mass, momentum, induction, and entropy equations are re- 
spectively 

^ + V.(pu) = 0, (1) 

Dv (VxB)xB ^ n , o . 

Pm-to VP + m (2) 

Vx(vxB), (3) 



dt 
D In Pp~^ 7 - 1 



[V-Q + pC], (4) 



Dt P 

where p is the mass density, v is the fluid velocity, B is the magnetic field, g is the gravita- 
tional acceleration, P is the gas pressure, 7 is the adiabatic index (5/3 for monotonic gas), 
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Q is the heat flux, and C is the radiative energy loss per unit mass of fluid, whose form we 
will leave unspecified. For thermal bremsstrahlung, a reasonable approximation is 

pC ~ 2 x 10~ 27 n 2 T 1/2 ergs cm" 3 s _1 , (5) 

where n e is the electron number density. D/Dt is the Lagrangian derivative, d/dt + v-V. 

To define the heat flux Q, let 6 be a unit vector in the direction of the magnetic field. 
Then (Balbus 2001): 

Q = - X b(b.V)T, (6) 

where T is the kinetic gas temperature and \ is the thermal conductivity (Spitzer 1962): 

X ~ 6 x 10~ 7 T 5/2 ergs cm" 1 s" 1 K _1 . (7) 

Finally, we follow Q08 and use the notation 

k = X T/P (8) 

for the thermal diffusion coefficient. 

2.1. Equilibrium Background 

We consider a gas stratified in the vertical z direction, with temperature profile T(z). 
The gas is not self-gravitating, so that g is a specified function of position. We assume a 
highly sub-thermal magnetic field. Thus, in equilibrium, the gas is in hydrostatic balance, 

dP 

^ = -«. (9) 

The magnetic field is uniform with x and z components B x and B z , (in this way defining the 
x axis), and unit vectors b x = B x /B, b z = B z /B. In equilibrium, there is a thermal balance 
between conductive heating and radiative losses, 

-v.g s « = ^ do) 

2.2. Local WKB Perturbations 

As in Q08, we consider plane wave disturbances of the form exp(a£ + ik • r) where the 
wavenumber k has Cartesian components (k x , k y , k z ), and r is the position vector. We differ 
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in notation from Q08 by using a, a formal growth rate, rather than uj, an angular frequency. 
This ensures that all coefficients in the final dispersion relation are real. We work in the 
WKB (At 3> 1) and Boussinesq limits (Q08). 

We next consider the linearized equations when perturbations Sp, Sv, etc., are applied 
to the equilibrium state. The heart of the problem is the entropy equation, so let us begin 
here. The linearized form of equation (j4j) is 



P 



where 



and 



rflnPp-T 

dz 
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V-SQ 
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= pC/P, 
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■ 
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(13) 



that is, the derivative of G with respect to T with P held constant. We have used the 
Boussinesq approximation in ignoring all terms proportional to 6P in equation ( ITT]) . In the 
process, we have implicitly regarded as a function of T and P (rather than the more 
customary but less convenient T and p dependence). The remaining linearized equations, 



_6pdP_ 
p 2 dz 



k-v 

V p 



= 0, 
B5B 



i(k ■ B)5B 
4np 



(14) 

(15) 
(16) 



4irp 
a5B = i(k • B)Sv, 

are, apart from notational convention, identical to Q08. The entire system of equations 
differs from Q08 only by the G term. The resulting dispersion relation is 



a 



7 



7 



TQ Tl p + C) (a 2 + (k-v A ) 2 ) 



where 



N 2 



ak\N 2 „ « dlnT 
+ h + C/C^,^— = 0, 



k 2 k 2 dz 



1 dPdlnPp-~< 



di n p(i-7)/7 T 



P7 dz 
C = 



dz 

7-1 

7 



g- 
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n{k-b)' 
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(19) 
(20) 
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This corresponds to equation (13) of Q08, except, as noted, for the single appearance of the 
radiative &t\p term. (A less general version of this result was also presented in Balbus & 
Reynolds [2008].) The dispersion characterizes the linear response of a magnetized, thermally 
conducting radiative dilute plasma to incompressible disturbances. 



2.3. Stability 

2.3.1. Recovery of the Conductive Field criterion 



Expanding the dispersion relation (jT7|) leads to 

a 3 + aicr 2 + a 2 cr + a 3 = 0, (21) 

where 

01= (^—^\Te Tl p + C, (22) 

a 2 = ^N 2 + (k ■ v A )\ (23) 

a 3 = CK -^~^~ + ( " Va) V (24) 

There are stable solutions to this dispersion relation if and only if the following three criteria 
are met: 

Oi > 0, a 3 > 0, a\a2 > a 3 . (25) 

This follows from a Routh-Hurwitz analysis, but can be seen more easily by inspection: the 
first two are in fact elementary, while the third follows from self-consistently demanding 
purely imaginary solutions to the cubic equation and then investigating their behavior for 
infinitesimal real parts. 

The physical interpretation of a\ > 0, or 

T6 T |p + (k-b) 2 K > 0, (26) 

is the magnetized conduction variation of the classical thermal instability criterion (Field 
1965). Only the component of k along the field lines enters into the conduction term. 



2.3.2. Recovery of the HBI and MTI 

We next consider the physical interpretation of a 3 > 0, or 

„„„ q d\n.T n9 . . 

c ¥~dT Jr{ - a) 1>0 " (27) 
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In essence, this is the HBI/MTI criterion of Q08, but further (de) stabilized when the flow 
is (un)stable by the isobaric Field criterion. This is a true instability if a 3 is negative, with 
a = —a 3 /a 2 in the limit of large a 2 > 0. 

Equation (F27I) shows that thermal instability and the HBI/MTI are intimately linked. 
To be definite, consider the behavior of the HBI. The discussion of Q08 explains how the 
distortion of field lines leads to conductive cooling of a downwardly displaced fluid element 
(and vice- versa for an upwardly displaced element). It is this cooling that causes the con- 
vection associated with the HBI. With radiative losses present, the cooling is enhanced on a 
downward displacement, and relative heating is present on an upward displacement. In fact, 
we may imagine now slowly turning on the magnetic field from dynamically weak to strongly 
dominant. Then, the exponentially growing instability transforms from the Q08 HBI to the 
classical (nonoscillatory) thermal instability. The role of the magnetic field in mediating this 
transition is crucial. 



2.3.3. Destabilization of wave modes by a positive thermal gradient 

We now return to the third criterion, a\a 2 > a 3 . With a% > 0, this criterion is a more 
stringent stability criterion than the first (oi > 0), and hence replaces it. 

When a 2 is large and positive [e.g. either N 2 or (k • va) 2 is dominant], the unstable 
roots depending on a± will be approximately: 

1 /9 

a = ±ia 2 + («3 — aia2)/2a2 (28) 

On the other hand, at large wavenumbers, we may have a\ and a% as the dominant terms. 
If a?, and a-i are both positive (or both negative), then the wavelike solutions will be 

a = ±i(a 3 /a 1 ) 1/2 + (o 3 - a 1 a 2 )/2a 2 1 (29) 

In either case above, the combination 03 — a\a 2 determines the stability of the mode. 

After a cancellation of the magnetic tension terms, the condition a\a 2 — a 3 > becomes 

k 2 \ 9 „,^ Q dlnT , . 

^-^ihr >0 - (3 °) 

Consider the limit 6 Z <1, which is HBI stable (03 > 0) for all but nearly axial wavenumbers, 
whose growth times then become very long. (We cannot take b z = exactly, since that 
would preclude a static radiative equilibrium state. For b z finite, equation (ITUl) shows that 
the equilibrium dT/dz scales as b~ x .) Then, /C = kj_, and our inequality reduces to 

7 — 1 r^^ , T 9 „ /, T 9 dlnT 



TQ T]P N 2 + C[N z -g——)>0. (31) 

7 V dz ' 



But 

and assuming that N 2 > 0, the inequality may be yet further reduced: 

C fdP x 2 



Te ^ + jpwA^) >0 - (33) 

Finally, substituting for C and N 2 and simplifying, our condition becomes 

TQ T \ P + n(k ■ bfTZ > 0, (34) 

where 1Z is the reduction factor 



K 



7 d\nT N ' 



7-ldlnP 



(35) 



which bears direct comparison with equation ( )26|) . Here, thermal conduction once again 
stabilizes radiative losses, but the HBI terms, when they are in a stable configuration rel- 
ative to the convective processes discussed by Q08, actively destabilize wave-like modes by 
reducing the supression of thermal conduction. In regions of sharp temperature gradients, 
the effective reduction factor for conductive stabilization can be large. Indeed, in the chosen 
limit b z <C 1, we have TZ ~ 0(b z ). Note that wavenumbers with vanishing k • b are unaf- 
fected by conduction, and have an effective reduction factor of zero. In our example, these 
are horizontal fluid displacements along the magnetic field lines. 



2.3.4- Destabilization of wave modes by a negative thermal gradient 

Consider next the case b z — 1, which would be HBI- unstable in the case of an increasing 
outward temperature profile. But let us now assume that the temperature decreases out- 
wards. This configuration is HBI stable. With b z = 1, if we restricted ourselves only to the 
first two stability criteria, we would conclude that this configuration is also MTI-stable. In 
fact, if the third stability criterion is imposed, this configuration is subject to an interesting 
and powerful overstability, driven by anisotropic thermal conduction, as we now show. 

With b z = 1 we have /C = — k\ and our third criterion inequality (|30|) becomes 

?—± TQ TlP N 2 + C (n 2 + g^-) > 0. (36) 

7 V dz J 
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Once again, the thermal conduction is affected by a "reduction factor," though here the 
reduction factor 1Z' actively destabilizes rather than merely supresses dissipative destabiliza- 
tion. The above inequality may be written 

TQ T \ P + «(fe • b) 2 W > 0, (37) 

where 

W = 1- ''-'— 1 (38) 



7 dlnT 
The term inside the square brackets must always be positive if N 2 > 0, but if 

7 - 1 d In P 

1< —7^ <2 < 39 > 

then 7Z' < and buoyant modes are overstable, even if there is no radiative loss term. 



2. 3. 5. Summary 

The MTI and HBI are evansecent instabilities present in dilute plasmas when anisotropic 
heat flux is included in the physics. The MTI is present when the thermal gradient decreases 
outward and the field lines are insulating in the equilibrium configuration. When the field 
lines open, the MTI is stabilized. The HBI is present when the thermal gradient increases 
outward and the field lines are open so that a heat flux is present in the equilibrium config- 
uration. The action of the HBI is to close the field lines, which stabilizes the system. 

We have found that the stable "end states" of these instabilities are subject to further 
overstabilities. In the case of the HBI, which is relevant for the coooling flow cores, a ther- 
mally unstable radiative loss function and closed fields lines together manifest as over stable 
buoyant oscillations. In the case of the MTI, a sufficiently steep (but classically convectively 
stable) outwardly decreasing thermal gradient produces overstable buoyant waves when the 
magnetic field lines are open and conducting heat. 

The overstabilities nominally depend on radiative losses, but their effect should be 
thought of as dynamical: these are classical g-waves that in principle could be driven to 
finite amplitudes on thermal time scales (either radiative or conductive). Whether they are 
best thought of a local WKB waves, global modes, or both is not yet clear, and awaits 
numerical investigation. 
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3. Discussion and Conclusions 

The implications of the Q08 finding that generic cluster (or elliptical galaxy) cooling 
flows are convectively unstable have yet to be grasped. A more complete linear theory 
is clearly a starting point. Here, we have generalized the linear theory of such systems to 
include the effects of both anisotropic thermal conduction and optically-thin radiative losses. 

To recap, strict stability requires three criteria to be satisfied. The first amounts to the 
classical Field criterion for thermal instability in the presence of anisotropic conduction, 

d = TQ T \ P + n(k-b) 2 > (Stability). (40) 

The second criterion gives the MTI or HBI stability conditions depending upon the orienta- 
tion of the magnetic field (via /C) and the temperature gradient, 

clc gan_ + , fe e x 2(li > (Stability). (41) 

k 2 dz 

The third criterion has not, to our knowledge, been recognized previously. For the two 
limiting cases considered in this work, it takes the form 

re T | P + K(k-bfTZ > (Stability ; b z « 0; dT/dr > 0), (42) 

TQ T \ P + K(k-b) 2 TZ' > (Stability ; b z w 1; dT/dr < 0), (43) 

where < TZ < 1 and — oo < 1Z' < 1. Even once the HBI (MTI) has been stabilized by the 
formation of horizontal (vertical) magnetic fields during their non-linear evolution, the third 
criterion can be violated in some range of wavenumbers leading to overstable (/-modes. 

The mechanism of the overstability differs depending upon the setting. The overstability 
that would naturally follow HBI-driven evolution (dT /dz > and essentially horizontal 
field) is driven by the presence of Field-unstable radiative cooling. On the other hand, the 
overstability that would naturally follow MTI-driven evolution (dT ' /dz < and essentially 
vertical field) is driven by the anisotropic heat flux, and is present even in the absence of 
radiative losses. 

Figure 1 indicates schematically the set of instabilities and overstabilities discussed in 
this paper. At a formal level, it is remarkable that the conductive heat flux can drive the 
atmosphere away from stability irrespective of the sign of the temperature gradient: only very 
shallow negative temperature gradients [such that dlnP/dlnT > 27/(7 — 1)] are completely 
stable, and even here the overstability can be reinstated if Field-unstable radiative losses are 
present. 
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Fig. 1. — Schematic map of the instabilities and overstabilities discussed in this work. 

The astrophysical implications of these overstabilities have yet to be determined, and 
we must await simulations in order to assess their global character as well as their non-linear 
behavioufl Here, we limit ourselves to some brief remarks and speculation. Probably tbe 
most important implication of these overstabilities for cooling core clusters is that, provided 
that radiative losses maintain a temperature gradient, the ICM may never be able to achieve a 
state that is stable to the HBI and the g-mode overstability. This statement holds true even 
if some unspecified heat source is fending off core collapse, provided that the net cooling 
function is still Field-unstable. The g-mode overstability grows on the cooling timescale 
and can be easily seeded by perturbations resulting from AGN activity, galaxy motions, 
sub-cluster mergers or indeed a preceeding phase of HBI-driven turbulence (Ruszkowski 
& Oh 2010). The g-mode overstability may still be relevant even when forced turbulence 
erases any obvious manifestations of the HBI. It is interesting to speculate that the "cold 



1 To date, published local studies of the HBI have excluded radiative cooling, and so are not be subject 
to this overstability. The overstabilities should, however, be present in published global studies of cooling 
cores (Bogdanovic et al. 2009; Parrish et al. 2009). But these simulated cores undergo a thermal runaway 
after several cooling times, and the overstabilities may be buried in the complex background dynamics, and 
very difficult to extract. 
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fronts" seen in the ICM cores of many relaxed clusters (Markevitch et al. 2000; Ascasibar & 
Markevitch 2006) may be related to the non-linear development of these radiatively-driven 
(/-mode overstabilities. 

At the purely formal level, understanding the development of thermal instability in a 
dilute plasma (Field 1965) has turned into a four decade struggle. And one speaks here 
of linear stability! It is a problem of some subtlety. With the advent of the HBI and our 
deepened understanding of the role of anisotropic conduction, we may at last have the true 
essence of the problem. There may, on the other hand, be more to come. 



Acknowledgements 

CSR thanks the Physics Department of the Ecole Normale Superieure de Paris for its 
hospitality and support of a one month visit in June 2010 during which some of this work 
was conducted. This work is supported by the NSF under grant AST 09-08212. 



REFERENCES 

Ascasibar, Y., & Markevitch, M. et al., 2006, ApJ, 650, 102 

Balbus, S.A. 2000, ApJ, 534, 420 

Balbus, S.A. 2001, ApJ, 562, 909 

Balbus, S.A. 2004, ApJ, 616, 857 

Balbus S.A., & Reynolds, C.S. 2008, ApJ, 681, L65 

Betschinger, E., & Meiksen, A. 1986, ApJ, 306, LI 

Bogdanovic, T., Reynolds, C.S., Balbus, S.A., Parrish, I.J., 2009, ApJ, 704, 211 

Braginskii, S.I. 1965, Rev. Plasma Phys., 1, 205 

Conroy, C, & Ostriker, J.P, ApJ, 681, 151 

Field, G.B. 1965, ApJ, 142, 531 

Markevitch, M. et al., 2000, ApJ, 541, 542 

Parrish, I.J., & Quataert, E. 2008, ApJ, 677, L9 



-13- 

Parrish, I.J., Quataert, E., Sharma, P., 2009, ApJ, 703, 96 

Parrish, I.J., Quataert, E., Sharma, P., 2010, ApJ, 712, L194 

Parrish, I.J., Stone, J.M. & Lemaster, N.. 2008, ApJ, 688, 905 

Peterson, J.R., & Fabian, A.C., 2006, Physics Reports, 427, 1 

Ruszkowski, M., Oh, S.P, 2010, ApJ, 713, 1332 

Quataert, E. 2008, ApJ, 673, 758 (Q08) 

Zakamska, N.L., & Narayan, R. 2003, ApJ, 582, 162 



This preprint was prepared with the AAS IATgX macros v5.2. 



